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ABSTRACT 

We present the first study on the mass functions of Jenkins et al (J01) and 
an estimate of their corresponding largest virialized dark halos in the Universe 
for a variety of dark-energy cosmological models with a running spectral index. 
Compared with the PL-CDM model, the RSI-CDM model can raise the mass 
abundance of dark halos for small mass halos at lower redshifts, but it is not ap- 
parent on scales of massive mass halos. Particularly, this discrepancy increases 
largely with the decrease of redshift, and the RSI-CDM model can suppress the 
mass abundance on any scale of halo masses at higher redshift. As for the largest 
mass of virialized halos, the spatially flat ACDM models give more massive mass 
of virialized objects than other models for both of PL-CDM and RSI-CDM power 
spectral indexs, and the RSI-CDM model can enhance the mass of largest virial- 
ized halos for all of models considered in this paper. So we probably distinguish 
the PL-CDM and RSI-CDM models by the largest virialized halos in the future 
survey of cluster of galaxies. 

Subject headings: cosmology:theory — dark matter — galaxies:halos — large-scale 
structure 



1. Introduction 

The central problem in modern cosmology is the formation of large scale structures 
in the universe. In the standard picture of hierarchical structure formation, dark matter 
dominates the universe, and a wide variety of observed structures, such as galaxies, groups 
and clusters of galaxies, have formed by the gravitational growth of Gaussian primordial 
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density fluctuations. Due to self-gravitational instability, the fluctuations of dark matter 
have collapsed and virialized into objects which are so-called 'dark matter halos' or 'dark 
halos'. The larger halos are generally considered to have formed via the merger of smaller 
ones collapsed first. The distribution of mass in the gravitationally collapsed structures, 
such as galaxies and groups (or clusters) of galaxies, which is usually called the mass or 
multiplicity function, has been determined by observation. 

As the observational data relevant to these issues improve, the need for accurate theo- 
retical predictions increases. By far the most widely used analytic formulae for halo mass 
functions are based on extensions of the theoretical framework first sketched by Press & 
Schechter (1974). The Press-Schechter (PS) model theory did not draw much attention until 
1988, when the first relative large N-Body simulation revealed a good agreement with it. 
The mystery of the 'fudge factor' of 2 in PS theory was solved by approaching the 'cloud- 
in-cloud' problem with a rigorous way(Peacock & Heavens 1990; Bond et al. 1991). The 
reliability of the PS formula has been tested using N-Body simulation by several authors, 
which turns out the PS formula indeed provides an overall satisfactory description of mass 
function for virialized objects. Unfortunately, none of these derivations is sufficiently rigor- 
ous such that the resulting formulae can be considered accurate beyond the regime where 
they have been tested against N-body simulations. Although the analytical framework of 
the PS model has been greatly refined and extended in recent years, in particular to allow 
predictions for the merger histories of dark matter halos (Bond et al. 1991), it is well known 
that the PS mass function, while qualitatively correct, disagrees in detail with the results of 
N-body simulations. Specifically, the PS formula overestimates the abundance of halos near 
the characteristic mass and underestimates the abundance in the high mass tail. In order to 
overcome this discrepancy, Jenkins et al. (2001) proposed an analytic mass function which 
gives a fit to their numerical multiplicity function. 

In particular, a power spectrum of primordial fluctuation, P p (k), should be assumed in 
advance in the calculation of mass function. Inflationary models predict a approximately 
scale-invariant power spectra for primordial density (scalar metric) fluctuation, P p (k) oc k n 
with index n = 1 (Guth & Pi 1982; Bardeen et al. 1983). The combination of the first- 
year Wilkinson Microwave Anisotropy Probe (WMAP) data with other finer scale cosmic 
background (CMB) experiments (Cosmic Background Imager [CBI], Arcminute Cosmol- 
ogy Bolometer Array Receiver [ACBAR]) and two observations of large-scale structure (the 
Anglo-Australian Telescope Two-Degree Field Galaxy Redshift Survey [2dFGRS] and Ly- 
man a forest) favour a ACDM cosmological model with a running index of the primordial 
power spectrum (RSI- ACDM), while the WMAP data alone still suggest a best-fit standard 
power-law ACDM model with the spectral index of n ss 1 (PL-ACDM) (Spergel et al. 2003; 
Peiris et al. 2003). However, there still exist the intriguing discrepancies between theoretical 
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predictions and observations on both the largest and smallest scales. While the emergence 
of a running spectral index may improve problems on small scales, there remain a possible 
discrepancy on the largest angular scales. It is particularly noted that the running spec- 
tral index model suppress significantly the power amplitude of fluctuations on small scales 
(Spergel et al. 2003; Yoshida et al. 2003). This imply a reduction of the amount of sub- 
structure within galactic halos (Zentner & Bullock 2002). Yoshida et al. (2003) studied 
early structure formation in a RSI-ACDM universe using high-resolution cosmological N- 
body/hydrodynamic simulations. They showed that the reduced small-scale power in the 
RSI-ACDM model causes a considerable delay in the formation epoch of low-mass minihalos 
(~ 1O 6 M ) compared with the PL-ACDM model, although early structure still forms hier- 
archically in the RSI-ACDM model. Thus the running index probably affect the abundance 
of dark halos formed in the evolution of the universe. 

Among the virialized structures, galaxy clusters are extremely useful to cosmology be- 
cause they may be in detail studied as individual objects, and especially are the largest 
virialized structure in the universe at present. The mass of a typical rich clusters is approx- 
imately 1O 15 /i _1 M , which is quite similar to the average mass within a sphere of 8/i _1 Mpc 
radius in the unperturbed universe. However, the theoretical estimate of the mass of the 
largest collapsed object in the RSI-ACDM cosmological framework has still not been pre- 
sented. Therefore, we will calculate the mass function of collapsed objects by J01 mass 
functions respectively and present the first calculation of the largest virialized object in the 
Universe in a RSI-ACDM model to explore the effect of running spectral index of primordial 
fluctuation on structure formation. 

The reminder of this paper is organized as follows. We describe mass function of dark 
halos in Section 2. The largest virialized dark halos in the universe are presented in Section 
3. The conclusion and discussion are given in Section 4. 



2. Mass Function of Dark Halos 

In the standard hierarchical theory of structure formation, the comoving number density 
of virialized dark halos per unit mass M at redshift z can be expressed as: n(M, z) = 
dN/dM = p f(M,z)/M where p is the mean mass density of the universe today and, 
instead of PS formula in this letter, the mass function /(M, z) takes the form of an empirical 
fit from high- resolution simulation (Jenkins et al. 2001) 

. 0.301 din a' 1 (M,z) . „ . ,oo Sn 

f(M,z) = — rfln M exp(-|ln O r- 1 (M,^) + 0.64| 3 - 88 ). (1) 
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Here a(M,z) = a(M)D(z) and D(z) = e(£l(z)) /e(£l m )(l + z) is the linear growth function 
of density perturbation (Carroll et al. 1992), in which e(x) = 2.5x/(l/70 + 209x/140 — 
x 2 /140 + a; 4 / 7 ) and Q(z) = O m (1 + z) 3 / E 2 (z) . The present variance of the fluctuations within 
a sphere containing a mass M can be expressed as <x 2 (M) = P(k)W 2 (krM)k 2 dk, 
where J¥(/crM) = 3[sin(/crM)/(fcrM) 3 — cos(/crM)/(^M) 2 ] is the Top-hat window function in 
Fourier space and r M = (3M / 'Airpo) 1 ^ . The power spectrum of CDM density fluctuations 
is P(k) = P p (k)T 2 (k) where the matter transfer function T(k) is given by Eisenstein & 
Hu (1999), and P p (k) is the primordial power spectrum of density fluctuation. The scale- 
invariant primordial power spectrum in the PL-ACDM model is given by P p (k) = Ak n " with 
index n s =l and that in the RSI-ACDM model is assumed to be P p (k) = P(k )(k/k ) n3( - k \ 
where the index n s {k) is a function of length scale 

The pivot scale /c =0.05 h Mpc -1 , n s (/c )=0.93, and dn s /d\n /c=-0.03 are the best-fit values to 
the combination data of the recent CMB experiments and two other large-scale structure ob- 
servations (Spergel et al. 2003). For both PL-ACDM and RSI-ACDM models, the amplitude 
of primordial power spectrum, A and P(k ), are normalized to cr 8 = o"(r M = 8/i~ 1 Mpc), which 
is the rms mass fluctuations when present universe is smoothed using a window function on 
a scale of 8/i _1 Mpc. In this section, we assume spatially flat ACDM models characterized 
by the matter density parameter Q m , vacuum energy density parameter Q\. For both PL- 
ACDM and RSI-ACDM models, we take cosmological parameters to be the new result from 
the WMAP: Hubble constant h = 0.71, Q m = 0.27, a 8 = 0.84 (Bennett et al. 2003; Spergel 
et al. 2003). 

The mass function of dark halos directly involve the calculation of primordial power of 
density fluctuation. In order to explore the difference between the two kinds of primordial 
power spectrum, we first calculate the mass function of dark matter halos in a wide range 
of redshift, which are plotted in Fig.(l). It is noted that there is a slight difference between 
the PL-CDM model and RSI-CDM model at lower redshifts. Compared with the PL-CDM 
model, the RSI-CDM model can raise the mass abundance of dark halos for small mass halos 
at lower redshifts, but it is not apparent on scales of massive mass halos. Particularly, this 
discrepancy increases largely with the decrease of redshift. Similar to the result (Yoshida 
et al. 2003), the RSI-ACDM model can suppress the mass abundance on any scale of halo 
masses at higher redshift. According to the hierarchical formation theory of structure, there 
is fewer higher mass halos at higher redshift and the higher mass halos are formed by the 
merger of lower mass haloes at the relative late stage. As pointed previously, the RSI model 
can suppress the power spectrum at small scale, so this just leads to a considerable delay in 
the formation of low mass haloes instead of high mass haloes. Therefore, compared with the 
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PL-CDM model, the mass function is uniformly lower for the RSI model at higher redshift 
of z—6. 



3. The Largest Virialized Dark Halos In the Universe 

Based on the theoretical expression above, we can easily get the total number N of the 
virialized objects with the mass larger than M 



N 



OO POO 



J M 



dn 
dM 



dV 
dM]^-dz, 

dz 



(3) 



where dV is the comoving volume element for the Friedmann-Robertson- Walker metric and 
^j- takes the form 

dz 



dV 

dz 



H o y/l+n k D*(l+z)*E(z) 
c 3 Dl{l+zf 



— = < 4n M 



E(z) 



for Q k = 



(4) 



i c 3 .D„ 2 (l+*) 2 cosh[^./] f Q, > n 



4 y/\+n k Dl{\+z)*E{z) 

where D a = dAH /c, dA is the angular diameter distance and / = J Q Z dz/E(z). 

It is obvious from the Eq. (3) that the total number N decrease with the increase of the 
mass M. Setting N = 1, we can finally obtain the largest mass Mm ax of virialized object 



dn .^(/V , 
dM]— dz. 



M M ax 



dM 



dz 



(5) 



In this section, we consider three cold matter (CDM) models, i.e. the standard CDM 
(SCDM), spatially flat ACDM models, and an open CDM (OCDM)for both PL and RSI 

Table 1. Cosmological Models Parameters 



Model Q m Q, 



OS 



SCDM 1 0.5 0.6 
LCDM 0.3 0.7 0.21 1.0 
OCDM 0.3 0. 0.25 1.0 
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Fig. 1. — Mass function at redshift z—0, 0.6, 1.27 and 6, respectively. The solid line is the 
mass function for running spectral index ACDM model, while the dashed one is that for 
power law ACDM model. 
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spectrum models. The cosmological models parameters are given in Table 1. Then we cal- 
culate the largest virial mass Mmax i n a variety of cosmological models for both PL and 
RSI power spectrum model, the results of which are demonstrated in Table 2. From Table 2 
we can see that the different cosmological models may yield the different result about virial 
mass for the largest virialized halos. The spatially flat ACDM models give more massive 
mass of virialized objects than other models for both of PL and RSI power spectral models. 
Therefore, it can distinguish different cosmological models by the largest mass of virialized 
halos. Due to the accumulative effect of the integration for volume(or redshift) over the 
whole space in the universe, the prediction for virial mass is slightly greater than the ob- 
served typical one. In addition, we also notice that the RSI-CDM model can enhance the 
mass of largest virialized halos for all of models considered here. 



Motivated by the new result on the index of primordial power spectrum from a combina- 
tion of WMAP data with other finer scale CMB experiments and other large-scale structure 
observations, we present the first study on the mass functions of J01 and their correspond- 
ing largest virialized dark halos in the Universe for a variety of dark-energy cosmological 
models with a running spectral index. It is well known that structures in the universe forms 
hierarchically in standard CDM models. The most massive structure form rather late in 
the universe. It is also noted that there is a slight difference between the mass abundance 
of PL-CDM and RSI-CDM model at lower redshifts. Compared with the PL-CDM model, 
the RSI-CDM model can raise the mass abundance of dark halos for small mass halos at 
lower redshifts, but it is not apparent on scales of massive mass halos. Particularly, this 
discrepancy increases largely with the decrease of redshift, and the RSI- ACDM model can 
suppress the mass abundance on any scale of halo masses at higher redshift. 

Table 2. Numerical Results for the largest virial mass Mm j 4x(10 15 /i~ 1 M q ) 
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Conclusions and Discussion 



Model PL model RSI model 



SCDM 
LCDM 
OCDM 



3.0 
6.2 
4.5 



3.3 
6.9 
5.05 
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As for the largest mass of virialized halos, the spatially flat ACDM models give more 
massive mass of virialized objects than other models for both of PL-CDM and RSI-CDM 
power spectral models. Therefore, it can distinguish different cosmological models by the 
largest mass of virialized halos for both of PL-CDM and RSI-CDM models. In addition, we 
also notice that the RSI-CDM model can enhance the mass of largest virialized halos for all 
of models considered here. So we probably distinguish the PL-CDM and RSI-CDM models 
by the largest virialized halos in the future survey of cluster of galaxies. Therefore, the 
obtained largest virialized object can be referred to as the complement to the observations 
of CMB, SN la and large scale structure in the future cosmological observation. 

Yoshida et al. (2003) found that although the hierarchical formation mechanism do not 
work well in RSI-ACDM model compared with that in PL-ACDM model and it also is not 
clear that the PS theory can be used in RSI-ACDM model, the mass function measured 
by high-resolution cosmological N-body/hydrodynamic simulations overall match the PS 
mass function for both RSI-ACDM and PL-ACDM model. In addition, because the running 
spectral index model predicts a significant lower power of density fluctuation on small scales 
than the standard PL-ACDM model(Spergel et al. 2003; Yoshida et al. 2003), it should also 
attract considerable attention in studies on strong lensing (Zhang et al. 2004; Chen 2003a,b, 
2004a,b; Zhang 2004) and weak lensing by large-scale structure(Ishak et al. 2004), especially 
on skewness(Pen et al. 2003; Zhang et al. 2003; Zhang & Pen 2005), which characterizes the 
non-Gaussian property of n field in the nonlinear regime. 

We are very grateful to the anonymous referee for valuable comments that greatly 
improved the paper. We also thank Ue-Li Pen for useful discussions. T.J.Zhang would like 
to thank CITA for its hospitality during his visit. This work was supported by the National 
Science Foundation of China (Grants No. 10473002 and 10273003), the Scientific Research 
Foundation for the Returned Overseas Chinese Scholars, State Education Ministry. 
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